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Abstract 

Using the backstepping design, we achieve exponential stabilization of the coupled Saint-Venant-Exner (SVE) PDE model of water 
dynamics in a sediment-filled canal with arbitrary values of canal bottom slope, friction, porosity, and water-sediment interaction under 
subcritical or supercritical fiow regime. The studied SVE model consists of two rightward convecting transport Partial Differential Equations 
(PDEs) and one leftward convecting transport PDE. A single boundary input control (with actuation located only at downstream) strategy 
is adopted. A full state feedback controller is firstly designed, which guarantees the exponential stability of the closed-loop control system. 
Then, an output feedback controller is designed based on the reconstruction of the distributed state with a backstepping observer. It also 
guarantees the exponential stability of the closed-loop control system. The fiow regime depends on the dimensionless Eroude number 
Fr, and both our controllers can deal with the subcritical {Fr < 1) and supercritical {Fr >1) fiow regime. They achieve the exponential 
stability results without any restrictive conditions in contrast to existing results. 
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1 Introduction 

Balance laws are the key point for modeling complex physi¬ 
cal systems that involve fluid mechanics, reactions, heat and 
mass transfer phenomena. In fluid mechanics, fundamen¬ 
tal balance equations expressing the conservation of certain 
quantities, such as the energy, the mass or the momentum in 
physical processes, lead to spatio-temporal differential equa¬ 
tions that express transport or diffusion phenomena. Such 
equations are the starting point for the design of various 
controllers that ensure the stability and operability of many 
engineering applications. Among those applications, we are 
interested in the stabilization of the hyperbolic SVE PDEs 
describing the flow and the bed evolutions in an open chan¬ 
nel [11,12]. The SVE model has attracted considerable at¬ 
tention over the past decades. Several theoretical and exper¬ 
imental studies have been proposed in the literature, consid¬ 
ering the flow and sediment characteristics of the water mo¬ 
tion. These studies also addressed the influence of the par¬ 
ticle size, shape and density. However, the control of such 


^ Corresponding author: A. Diagne. 

Email addresses: ababacar. diagne@it .uu.se (Ababacar 
Diagne), mdiagne@ucsd.edu (Mamadou Diagne), 
sht015@ucsd.edu (Shuxia Tang), krstic@ucsd.edu (Miroslav 
Krstic). 


systems, modeled by nonlinear hyperbolic PDEs, is left out 
in most of the studies. To the best of the authors’ knowl¬ 
edge, there are only a few results on the stabilization of SVE 
model in the existing literatures. 

Several strategies have been developed to control the flow 
dynamics in classical irrigation canals without the sediment 
layer during the last decades. We refer the reader to [22] 
in which the classiflcation of control problems and related 
methodologies is fairly addressed. Basically, the main pur¬ 
pose is the regulation of the water level at a desired height by 
adjusting the opening of the gates at the ends of the channel 
as boundary actuators. Eor instance, the synthesis of LQ con¬ 
troller can be found in [1,21,31], whereas [20,23] has stud¬ 
ied an Hoc control approach. Through semigroup approach, 
[32] proposed an integral output feedback controller using a 
linear PDE model around a steady state. Lyapunov analysis 
is investigated [5,24,14] and multi-models approach with a 
stability analysis based on Linear Matrix Inequality {LMI) 
is presented in [13,25]. 

Very recently, [29] proposed a singular perturbation ap¬ 
proach for the synthesis of boundary control for hyperbolic 
systems. The effectiveness of the controller is illustrated us¬ 
ing the linearized SVE model. In [11], explicit boundary 
dissipative conditions are given for the exponential stability 
in L^-norm of one dimensional linear hyperbolic systems 
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of balance laws. The proposed Lyapunov approach is ap¬ 
plied to the linearized SVE equations with successful results. 
However, on-line measurements of the water levels at both 
ends of the spatial domain, namely, at upstream {x = 0) and 
downstream {x = L) are assumed to be available. Later on, a 
priori estimation technique and the Faedo-Galerkin method 
are proposed in [12] for the design of a linear feedback con¬ 
trol law that requires only downstream measurements. Such 
an approach have been vividly presented in [10,15] as well. 


controller that guarantees the exponential stability property, 
requires a measurement at the two boundary of the channel 
while in the present work we just use the downstream gate 
for actuation without need of a sensor there. Better, with the 
backstepping method, we achieve the exponential stability 
around the origin without imposing any conditions on the 
matrix from the source term of the system under study. Con¬ 
versely in [11], that matrix is required to satisfy a restrictive 
condition to be marginally diagonally stable. 


The stabilization problems for hyperbolic systems have been 
widely studied in the literature. The first approach relies on 
careful analysis of the classical solutions along the charac¬ 
teristics. We refer the readers to [16] in the case of second- 
order system of conservation laws and to [19] for more 
general situations as ^th-order systems. Another approach 
based on the Lyapunov techniques is introduced by [6] and 
improved in [5] where a strict Lyapunov function in terms 
of Riemann invariants is constructed and its time deriva¬ 
tive can be made negative definite by choosing properly the 
boundary conditions. The aforementioned Lyapunov func¬ 
tion is very useful to analyze nonlinear hyperbolic systems 
of conservation laws because of its robustness. We refer 
to [3,4,5,7,12,30,33], in which several applications are ana¬ 
lyzed based on this tool. 

Recently, the backstepping method was introduced for the 
feedback stabilization of various classes of PDEs [9,18]. 
The key idea of this approach is the construction of suitable 
Volterra integral transformations that map the original sys¬ 
tem into a so-called “target system”, which is exponentially 
stable. The kernel functions of the transformations are re¬ 
quired to satisfy some PDEs, and the solutions can be then 
used as gains of the controllers. The invertibility of the trans¬ 
formations ensures the exponential stability of the closed- 
loop control systems. One can refer to [8,18,26,2] for further 
applications of this technique to other classes of systems in¬ 
cluding nonlinear PDEs. 

In the present work, we achieve an exponential stabilization 
of the coupled Saint-Venant-Exner (SVE) model of water 
dynamics in a sediment-filled canal with arbitrary values of 
canal bottom slope, friction, porosity, and water-sediment 
interaction [11]. The studied SVE model consists of two 
rightward and one leftward convecting transport PDEs [9]. 
First, a single boundary input control strategy (with actua¬ 
tion located only at downstream) is adopted and a full state 
feedback backstepping controller is designed to ensure the 
exponential stabilization of the system. Next, we employ a 
sole sensor at the upstream {x = 0) to derive an output feed¬ 
back controller based on the reconstruction of the distributed 
state with an exponentially convergent Luenberger observer. 
The fiow characteristics depend on the dimensionless Froude 
number, namely, Fr, and the proposed controllers operate 
under the subcritical (Fr < 1) and the supercritical (Fr >1) 
fiow regime. Moreover, the exponential stability results are 
obtained without imposing any restrictive conditions on the 
controller gains, which is in contrast to [11]. In this context, 
this paper is an extension of [11] in which the design of the 


This paper is organized as follows. In the next Section, the 
nonlinear SVE model is formulated based on its physical 
description, and a linearized version around a steady state is 
presented. Section 3 is dedicated to the backstepping trans¬ 
formation between the linearized model and a suitable ex¬ 
ponentially stable target system. Then, with the solutions 
to the gain kernel PDEs of the Volterra transformation, a 
full state controller is computed. An exponentially conver¬ 
gent backstepping observer is designed in 4. Based on the 
observer, which reconstructs the full state from the output 
measurement, an output feedback controller is constructed 
in section 5 and an exponential stability result is also es¬ 
tablished. Numerical simulations are provided for both the 
subcritical and supercritical fiow regime in Section 6, with 
a detailed discussion on the numerical computation of the 
controller gain. Finally in Section 7, a conclusion is pre¬ 
sented and some perspectives are discussed. 


2 The Saint-Venant-Exner model 

We consider a pool of a prismatic sloping open channel 
with a rectangular cross-section, a unit width and a moving 
bathymetry (because of sediment transportation). The state 
variables of the model are: the water depth H{t,x), the water 
velocity V{t,x) and the bathymetry B(t^x) which is the depth 
of the sediment layer above the channel bottom as depicted 
on Figure 1. The dynamics of the system are described by 
the coupling of Saint-Venant and Exner equations (see e.g. 
[17]): 


dH dH dV 

—\-V ^—\-H —— — 0 

ot ox ox 

dV dV dH dB 

ot ox 


(la) 

(lb) 

(l c) 


In these equations, g is the gravity constant, Sh is the bot¬ 
tom slope of the channel, C/ is a friction coefficient and a is 
a parameter that encompasses the porosity and viscosity ef¬ 
fects on the sediment dynamics. The coefficient a expresses 
as (cf [17]) 


where pg is the porosity parameter and Ag is the coefficient 
to control the interaction between the bed and the water fiow. 
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Figure 1. A sketch of the channel. 


2. 1 Steady-state and Linearization 


/ 0 0 o\ 

y* 

H* 

y 0 0 0/ 

The dimensionless Fronde number is defined as 


B{W*) = 


0 

y*2 

0 


C/77^ -2C/.^ 0 


Fr = 


y* 




Exact, but rather complicated expressions of the eigenval¬ 
ues of A{W*) can be obtained by using the Cardano-Vieta 
method, see [17]. Once the eigenvalues A/ of the matrix 
A{W*) are obtained, the corresponding left eigenvectors can 
be computed as 


A steady-state is a constant state (if*, V*, 5*)^ which sat¬ 
isfies the relation 


gSbH* = CfV*^. 


In order to linearize the model, we define the deviation of the 
state {H{t^x)^V{t^x), B{t,x)) with respect to the steady- 
state: 


(V* 


Lk = 


{Xk — Xi){Xk - Xj) 

e {1,2,3}. 




H*Xk 

gH* / 

(4) 


We multiply (3) by in order to rewrite the model in terms 
of the characteristic coordinates 1/4 (k = 1,2,3). Then we 
obtain 


^ h{x,t) ^ 
u{x^t) 
\^b{x,t)) 


V{x,t)-V* 
\B{x,t)-B* ) 


^ + A, ^ W, for /: = 1,2,3, (5) 

where 


Then the linearized system of the SVE model (1) around a 
steady-state is 


dh ,,^dh ^,^du 

^-l-y*^-l-if*^ =0 
dt dx dx 

dh db 


y* 


du ^j^du dh db V*^ . 


db 


^+aV*^^ = 0 . 


du 


dt 


dx 


( 2 a) 


( 2 b) 


( 2 c) 


^k = 


1 

{Xk — ^j) 


[{(V*-Xi){V*-Xj) 


+gH*)h + H*XkU + gH*b]. 


( 6 ) 


For the sake of simplicity, we introduce the following nota¬ 
tion rt: 

= r 

' fH*{X,-Xi){Xk-Xjy 

Some computations yields the following writing for equation 
(5): 


2.2 Characteristic (Riemann) coordinates 

In the matrix form, the linearized model (2) can be written 
as 

dW dW 

^+A(iy*)^=5(W*)W, (3) 

dt dx 

where 



(h\ 


(v* H* o\ 

w = 

u 

, A{W*) = 





0 aV*'^ 0 j 


f+A.f+ E(2 ;i,-3v>A = o. 

fork =1,2,3, (7) 

where the characteristic coordinates are now defined as 

^k = ( 8 ) 

rk 

Erom (7), the linearized model (5) in characteristic form can 
be written as 

f+ 4 -"«=»■ <’> 
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where 


then the system (13) is transformed into 


| = A = diag{X,MM). 


M 


^ a\ a2 as ^ 
ai a2 as 
\ ai a2 as ) 


with 

ak={3V*-2Xk)rk. 

From [17], the three eigenvalues of the matrix A are such 
that for a subcritical flow regime {Fr < 1), 


Ai < 0 < A 2 < As ( 10 ) 


and for a supercritical one {Fr >1), 

A 2 < 0 < Ai < A 3 (11) 


with Ai and A 3 being the characteristic velocities of the water 
flow and A 2 being the characteristic velocity of the sediment 
motion. Obviously, the sediment motion is much slower than 
the water flow. 


2.3 Change of notations 


Hereafter, we consider the case where the flow regime is sub- 
critical and adopt the following notations: v{t^x) = ^\{t^x), 
u\{t^x) = ^ 2 {t^x),U 2 {t^x) = <^ 3 (^v) and coefficients (charac¬ 
teristic velocities) Ai = —gi, Ji = A 2 and 72 = A 3 . We intro¬ 
duce also the vector n = (wi, ^ 2 )^^, the coefficients T7j = Oj+i 
for 7 = 1,2 and the matrix 


( a2 as 
\a2 as 


( 12 ) 


With the new variables, the set of equation (9) writes as: 

dtU\ FY\dxU\ = Oi\u\ F O 12 U 2 F aiv (13a) 

dtU2 + Yi^xUi = 02\u\ F O 22 U 2 F aiv (13b) 

dtv - jidjcV = riiuiF ri 2 U 2 + aiv. (13c) 


Introduce the variable 


w{t^x) = v(r,v)exp 



dtU\FY\dxU\ =G\\U\F0\2U2 

( a\ 

+ ai exp 1 ^ 


ai 

ai exp I —V 1w 


dtU2 + Y2dxU2 =G2\U\ F O 22 U 2 

Oi 

F 

dfW - jldjcW =r]i exp 

f ai \ 

-7772 exp ^ 2 - 


(14a) 


(14b) 


(14c) 


We rewrite this system as: 


dtui F Yi^xUi = G\iui + G\ 2 U 2 + a{x)w (15a) 

dtU 2 F Y 2 dxU 2 = G 21 U 1 + G 22 U 2 + a{x)w (15b) 

dfW — jldxW = B\{x)u\ F 02{x)u2 (15c) 


with a{x) = aiexp and 0j(x) = Oj+iexp for 

7 = 1 , 2 . 


To close the writing of the system (15), we enclose to it the 
following boundary and initial conditions 


uft^O) = qiw{t,0) fori =1,2, (16a) 

w{t, 1) = piui ( 7 , l)Fp 2 U 2 {t, l)FU ( 7 ), (16b) 

w(0,v) = w^(x), 7/i(0,x) = u^{x) for i = 1,2. (16c) 


Remark 1 Let us mention that in the case where the flow 
regime is supercritical, the following changes of variable 
will be considered (instead of the previous one) v{t^x) = 
^2{tF)y ui{t^x) = ^i{t,x), U 2 {t^x) = ^s{iF) and coefficients 
A 2 = —p, Yi = 2^1 72 = 2 . 3 . 


Ul(/,j) 



I-h 


x = 0 


1 


Figure 2. Schematic steep of the hyperbolic system. The internal 
coupling between the states in the system and boundary conditions 
are depicted. 
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As in [9], ui, U 2 and w are the distributed states and U(t) is 
the control input as shown in Figure 2. The measured output 
is given by: w(^,0) =y{t). 

3 Full State Controller Design 

3.1 Backstepping transformation and target system 
Consider the following backstepping transformation 





I-1 

;c = 0 1 


\l/i{t,x) = uft^x) for i = 1,2 (17) 

X{t,x)=w{t,x)- [ ki{x,^)ui{t,^)d^ 

Jo 

- [ k2{x,^)u2{t,^)d^- [ k3{x,^)w{t,^)d^. (18) 

Jo Jo 

We now seek a sufficient condition on the functions ki such 
that the transformation (17)-(18) maps the system (15)-(16) 
to the target system 

dtWi+YidxVi =(^nVi+(^i2W2 + a{x)x 

Jo 

+ r Ci2{x,^)W2{t,^)d^ 

Jo 

+ [ Ki{x,^)x{t,^)d^ (19a) 

Jo 

d, W2 + y2dxW2 =<721 + C722 V^2 + «WZ 

+ r C2i{x,^)Wi{t,^)d^ 

Jo 

+ r C22{x,^)W2{t,^)d^ 

Jo 

+ [ K2{x,^)x{t,^)d^ (19b) 

Jo 

dtX-lxdxX=0 (19c) 


Figure 3. Representation of the target system. 


A sufficient condition for the transformation (17)-(18) to 
map the original system (15) into the target system (19) is 
that the kernels ki satisfy the following system of first order 
hyperbolic PDFs: 

fldxki (x, ^) - 7 i ki (x, ^) 

= <7iiA:i(x,(^) + <72iA:2(x,(^) + ei((^)A: 3 (x,(*) (21a) 

fldxk2{x, ^) - 72 ^^ k2{x, ^) 

= <712^:1 (x, (^) + <722^2+ ^2('^)^3(-^) <?) (21b) 

ixdxki {x,^)+nd^ kj, (x, ^) 

= a{^)ki{x,^) + a{^)k 2 {x,^) (21c) 

with the following boundary conditions: 

k\{x,x) = k 2 {x,x) = (22a) 

n+M 72 + ^ 

pk 3 {x, 0 ) = qiYiki{x,0) ^q 2 Y 2 k 2 {x, 0 ). (22b) 

The existence, uniqueness and continuity of the solutions to 
the system (21) with boundary conditions (22) are assessed 
by Theorem 5.3 in [9]. 

Besides, plugging (17)-(18) into (19) and using (15)-(16) 
yields for / = 1, 2, 


with the following boundary conditions: 

V/i(r,0) = qiX{t,0) for / = 1, 2 and x{t, 1) = 0. (20) 


This dynamic is schematically represented on Figure . In the 
system (19), Cij(-) and K‘i(-) are functions to be determined 
on the triangular domain 


T= |(x,0 gM^|0<^ <x< i|. 

The system (19)-(20) is designed as a copy of the original 
dynamics with the coupling term in (15c) removed. As will 
be shown later, the new terms in (19a) and (19b) are neces¬ 
sary for the design but they will not affect the stability. 


0 = r\ai^)ki{x,^)-cnix,^) 

Jo 

+ / Ki{x,s)ki{s,^)ds ui{^)d^ 

J 

+ rfa(^)fc 2 (x,^)-c, 2 (x ,0 

Jo 

+ / K2{x,s)k2{s,^)ds U2{^)d^ 

+ rfa(^)fc3(x,^)-ifi(x,0 

Jo 

+ / Ki{x,s)k 3 {s,^)ds w{^)d^. 

J C d 


(23) 


(24) 


The coefficients iq can be chosen to satisfy the following 
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integral equation for / = 1, 2 


where for / = 1, 2, 


= a{x)k3{x,^)+ f Ki{x,s)k3{s,^)ds, (25) 

and the coefficients Cy can be chosen such that 

Cij {x,^) = a{x)kj{x,^) + Ki{x, s)kj(s,^) ds 

for/, 7 = 1,2 (26) 

under the fact that the ki exist and are sufficiently smooth. 


3.2 Inverse transformation and control law 


h{x,^)=ki{t,^)+ f ki{x,^)k{^,s)ds. (32) 

Thus, the control law U{t) can be obtained by plugging the 
transformation (18) into (15). Readily, xi^A) =0 implies 
that 


U{l)=-plUi{t,l)-p2U2{l,l) 

+ f \kiV,^)ui{x,^) 

Jo L 

+ k2{Wu2{x,^) + k2{l,^)w{l,^)\d^. 


(33) 


To ensure that the target system and the closed-loop system 
have equivalent stability properties, the transformation (17)- 
(18) has to be invertible. Since Xj/i = Ui, for i = 1,2, the 
transformation (18) can be rewritten as 


The ki in the integral term designate the kernel functions 
and satisfy the system ( 21 )-( 22 ). 


X{t,x)+ [ ki{x,^)\iri{t,^)d^ 

Jo 

+ rk2{x,^)W2{t,^)d^ 

Jo 

= w{t,x)- [ k 3 {x,^)w{t,^)d^. (27) 

Jo 

Let us define 

^{t,x)=x{t,x)+ [ ki{x,^)xi/i{t,^)d^ 

Jo 

+ r k 2 {x,^)yif 2 {t,^)d^ (28) 

Jo 


3.3 Stability of the target system and the closed-loop con¬ 
trol system 

We first prove exponential stability of the target system (19)- 

( 20 ). 

Lemma 1 For any given initial condition X^V ^ 

([0,1])) and under the assumption that Cy, /q G "^(T), 
the equilibrium (i/Ai, V/ 2 , = ( 0 , 0 , 0 )^ of the target 

system (19)-(20) is -exponentially stable. 


Since k^ is continuous by Theorem 5.3 in [9], there exists a 
unique continuous inverse kernel I 3 defined on T, such that 

w{t,x)=r{t,x)F ri 3 {x,^)r{t,^)d^, (29) 

Jo 

which yields the following inverse transformation Since 
Xj/i = Ui, for i = 1 , 2 , we could get the following relation 
from the first two equalities of (15) and (19): 

a{x)w = a{x)x+ f cn{x,^)Wi{t,^)d^ 

Jo 

+ [ ci 2 {x,^)W 2 {t,^)d^+ [ Ki{x,^)x{t,^)d^. (30) 

Jo Jo 

Thus, we could write the following 

w{t,x) = x{t,x)+ [ h{x,^)yfi{t,^)d^ 

Jo 

+ f l 2 {x,^)^ 2 {t,^)d^+ rh{x,^)x{t,^)d^, (31) 

Jo Jo 


PROOF. The stability proof is based on the time differen¬ 
tiation of the following Lyapunov function: 

V,it)=ta,e-^^f^ + ^)dx 
Jo V yi 72 J 

+ / f 7 -X^{t,x)dx, (34) 

Jo p 

where ai and 5i are strictly positive parameters to be deter¬ 
mined. 

Differentiating this function with respect to time, we get: 

VM=2 

Jo V 71 72 J 

+ 2 [ Fl^^dtXdx. (35) 

JO p 
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By taking into account the target system (19)-(20) and inte¬ 
grating by parts, we have 


Vi{t) = 


— a\e 

+ (1 ^ 

- / dx+ f 

Jo Jo 


x'P^(r,x)(-5i/2 + 2r,„,<7)»P(f,x) dx 

+ 2 [ aie~^^^'¥'^{t,x)ri„ya{x)xit,x) dx 

Jo 

+ 2 j j\^{t,x)ri^, 

X {C{x,^)'¥{t,^) +K{x,^)x{t,^ j) d^ dx, 


where the matrix a is defined in (12), the vectors 
a(v), K{x,^) and the matrices F/^v, C(x,^) are given by 


^{t,x) = 


K{x,^) = 
C(x,^) = 


M’’^A,aix)=H] 

\W 2 {t,x)J \a{x)) 


WxA))' ‘ 



/cn(x,^) Ci 2 (x,^) 

\c2l{x,^) C22{X,^) 


(36) 

(37) 

(38) 


Assume that for M > 0 and e > 0, we have 

||a||, |la(x)||,||C(x,^)||, \\K{x,m<M, 

Vxe[0,i],^e[0,x], (39) 

)^-(x)>e,V/=l,2,Vxe [0,1], (40) 


where the matrix/vector norms || • || are compatible with 
the other corresponding matrix/vector norms. Hence, using 
Young’s inequalities the following relations are derived 


2 f\ie-^^^'¥^{t,x)rinya'¥{t,x) dx 

Jo 

<2— f oie-^'^'P^(f,x)'P(f,x) t/x (41) 
£ Jo 


and 


2 r'I'^(r,x)ri„,C(x,0'P(f,^) dx 

Jo Jo 

M rx 

£ Jo Jo 

£ Jo 

+ ^ 1^1 (e-^>'*^-e-^>)'P^(f,x)'P(f,x)^/x 

<ai f (—x+'I'^(f,x)'P(f,x)Jx (43) 

Jo \£ oiS J 


2 [ aie [ ^*'^{t,x)rinyK{x,^)x{t,^) d^ dx 

Jo Jo 

<ail^ {'¥^{t,x)ri„yK{x,^) 

xK'^{x,^)rinv^{t,x)^X^{t:^))d^dx 

<( 2 i(^—^ J {t,x)^^{t,x)dx 

-^ai^ [ e~^^''x^{h^)dx. (44) 

Oi Jo 

Thus, using the boundary conditions (20), we obtain the 
following inequality 

- (^-ai{\ +dx 

-ax [\-^^’‘'¥'^{t,x)Px{x)'¥{t,x) dx, (45) 

Jo 

where 




2 r iny ^. 



h 


(46) 


2 [ aie (t,x)rinvCC(x)x(t,x) dx 

Jo 

< j( ('P^(?,x)r,>,va(x) 

xa^(x)r 

inv ^{t,x)^X^{t^^))dx 


I ^ {^^x)^{t,x)dx 



First, we choose the tuning parameter 5i > 0 sufficiently 
large so that the matrix P{x),x G [0,1] is positive definite. 
Then, by choosing 


0 <ai < min < 


5i 


2 2’5i + 1 

1=1 


(47) 


we could derive exponential stability of the target system. 
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with the boundary conditions 


Then, from the continuity and invertibility of the backstep- 
ping transformation (17)-(18), we could derive equivalence 
between the original system (15) (with the boundary and ini¬ 
tial conditions (16) and the control law (33)) and the target 
system (19)-(20). Thus, the following theorem is proved. 

Theorem 1 Consider the system (15) with the boundary and 
initial conditions (16) and the control law (33). Then under 
the assumptions that the initial data are in (^^([0,1]))^, 
the origin is exponentially stable in the sense. 

4 Backstepping Observer Design 

The feedback controller (33) requires a full state measure¬ 
ment across the spatial domain. In this section we are in¬ 
terested in the design of a boundary state observer for es¬ 
timation of the distributed states of the system (15)-(16) 
over the whole spatial domain using the measured output 
w(^,0) =y{t). The observer 


dfUi 

+ Jl^xUl 

= GuUi 

+ Cri 2 W 2 + oc{x)w 




-Pi{x) 

b(0 

-w(f, 0 )] 

(48a) 

dtU2 

+ yidxU2 

= G 21 UI 

H“ (722^2 C)C(x)w 




-Plix) 

b(0 

-w(f, 0 )] 

(48b) 

dtW 

- ndxW = 

- dl(x)u\ 

+ 02(v)w2 




-P3{x) 

b(0 

-w(f, 0 )], 

(48c) 


where (wi, ^ 2 , w)^ is the estimated state vector, consists of 
a copy of the plant plus an output injection and mimics the 
well-known finite dimensional observer format. The func¬ 
tions 0 j(v) = Oj+i for 7 = 1,2 and a{x) are the ones defined 
for the transformed system (15). The following boundary 
conditions have to be considered: 

Ui{t,0) = qiy{t) for i = 1,2 (49a) 

w{t, 1) = piui{t, l)^p 2 U 2 {t, l)-^U{t). (49b) 

Our objective is to find p\{x), P 2 {x) and P 3 {x) such that 
the estimated state vector (w, wi, U 2 ) converges to the real 
state vector (w, t/i, U 2 ) in finite time. Defining 

(^W ill ^2 ^ = (^W — W Ui —Ui U 2 — U 2 ^ (50) 

as the error variable vector, we obtain the following error 
system 

dfW — pdjcW = 6 i{x)ui + 62 {x)u 2 + P 3 {x)w{t, 0 ) (51a) 

dtUi^Yi^xhi =Gnui-\-(yi2U2-\-a(x)w 

-Cpi{x)w{t,0) (51b) 

dtiil + 72 dxU 2 =G 2 \U\ + G 22 U 2 + cx{x)w 

5 -p 2 {x)w{t,C) (51c) 


w(r,l) =piwi(r,l)-fp 2 M 2 (r,l), (52a) 

Wi(^,0)=0 fori =1,2. (52b) 


4.1 Backstepping transformation and the target error sys¬ 
tem 

Similarly to the controller design, we use the following in¬ 
vertible backstepping transformation 

iift.x) = ftft.x)^ [ mi{x,^)^{t,^)d^ (53a) 

Jo 

for i = 1,2 

w{t,x)=^{t,x)+ [ m 3 {x,^)^{t,^)d^, (53b) 

Jo 

where the kernels mf’) for /= 1, 2, 3 are defined in the 
triangular domain T to map the error system (51)-(52) into 
the following exponentially stable target system 

3~7ldx^l = Giifti + ( 7 i 2^2 

+ [ gnix,^)7ii{t,^)d^+ [ gi 2 {x,^)n 2 {t,^)d^ (54a) 
Jo Jo 

+ 72 <^ jc ^2 = ^21^1 + (522^2 

+ / g2l{x,^)7ti(t,^)d^+ [ g22{x,^)7t2{t,^)d^ (54b) 

Jo Jo 

dt^ - = 01 {x)7ti + 02(x)7r2 

-f / hi{x,^)7ti(t,^)d^^ [ h2{x,^)7t2{t,^)d^ (54c) 

Jo Jo 

with the boundary conditions 

iri(^,0) = 0 for / = 1, 2 (55a) 

0(r,l) =pi7ri(r,l)+p27r2(r,l). (55b) 

Here the functions gtj and hi have to be determined on 
the triangular domain 3^. As previously, we are attempting 
to find some sufficient condition for the kernels to match 
the target system. Differentiating the transformations (53) in 
time and space and substituting the results into (51) with the 
help of (54), the following PDEs are derived for the kernels 

Yidxmi — pd^mi = a\\m\ -f (7i2^2 + oc{x)m 3 , (56a) 
72 dxm2 - pd^m2 = (721^1 + ^ 22^2 + a(v)m3, (56b) 
pd^m^ -f pd^m^ = — 0i(v)mi — 02(v)m2. (56c) 
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To close the writing of the above system, the following 
boundary conditions are imposed: 


m\(x,x) = - a(x) 

n+M 

m 2 {x,x) = 


(57a) 

(57b) 

(57c) 


where 


n{x,^) =-mi{x,^) - f mi{x,s)r-i{s,^)ds. 


43 Stability of the target error system and convergence of 
the designed observer 


The observer gains are defined by 

Pi{x) = iimfx^O) for / = 1, 2, 3, (58) 

and the integral coupling coefficients are defined by the fol¬ 
lowing equations: 


hiix,^) = -e{^)m3ix,^) 


— j m3(x,s)/!,(s,i^)Pi', for i = 1, 2, 

(59a) 

Si,j{x,^) = -QA)m{x,^) 


- 1 mi{x,s)hj{s,^)ds, for {i, j} = 1, 2. 

(59b) 


4.2 Inverse Transformation 

The continuity of the kernel m3 in the transformation (53b) 
guarantees the existence of a unique continuous inverse ker¬ 
nel r3 in the transformation 

^{t,x)=w{t,x)+ f r-i{x,^)w{t,^)d^ (60) 

Jo 

define on T, and 


The observer is exponentially convergent to the original sys¬ 
tem. We first prove exponential stability of the target error 
system (54). 

Lemma 2 Under the assumptions \\f^, y/^, ^ 1]) 

and gij, hiG ^ (T), the system (54) with boundary conditions 
(55) and given initial condition exponentially 

stable in the sense. 

PROOF. The stability proof is based on the time differen¬ 
tiation of the following Lyapunov function 

h V n B / 

pi ^ 

+ / -0^(^,v)Jv, (63) 

Jo p 

where a 2 and 82 are strictly positive parameters to be de¬ 
termined. Differentiating this function with respect to time, 
we get: 

Jo V B » / 

A gS 2 X 

+ 2 / - ^dt^dx. (64) 

Jo p 


r 3 {x,^) = -m 3 {x,^)~ f m 3 {x,s)r 2 {s,^)ds. (61) 

Substituting (61) into (53 a), we obtain 

^{t,x) =Ui{t,x) - [ mi{x,^)w(t,^)d^ 

Jo 

- [ ( mi{x,^)rj.{^,s)w{t,s)dsd^ 

Jo Jo 

=Ui{t,x)- mi{x,^)w{t,^)d^ 

Jo 

~ lo dsd^, 


for / = 1, 2, 


and hence for / = 1,2, 


7ii{t,x) = Ui{t,x) + n{x,^)w{t,^)d^ (62) 


Taking into account of the target system (54) and integrating 
by parts, we rewrite (64) as 

yi{t)= \-a2e~^''{7tl{t,x)^7t2{t,x)) 




-52 


Jo 

+ 2/ a 2 e~^^Il^(t^x)TinyG Il(t^x)dx 
Jo 

— §2 / a 2 e~^^Il^(t^x) n(r,v) dx 

Jo 

+ 2 / f\ 2 e-^^lfit,x)rin, 

Jo Jo 

xG(v,(^) n{t,^)d^ dx 

A ^52X ^ 

4-2 / - ^(t^x)6(x)Il(t^x)dx 

Jo p 

A ^52X ^ nx 

+ 2 / - ^{t,x) / hix,^)Ilit,^)d^ dx, 

Jo p Jo 


dx 
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where the matrices O’ are defined by (38) and (12), the 
vector n, 0, A and the matrix G are given by 


no,,) = *■<'•'!), 

(65) 

\n 2 {t,x)) 



(66) 

\ g 2 l{x,^) g 22 {x,^)J 


0 {x) = (cc0i(x)) 62 {x)) , 

(67) 

h{x,^) = (yhi{x,^) h 2 {x,^) ) . 

(68) 


Assume that for M > 0, we have 


and finally 

/*! _ fx 

2 / -/ A(x,(^)n(r,(^) dx 

Jo 11 Jo 

pi ^d2X px y 

Jo /i Jo V 

+0^(r,x)) d^dx 

< J {t,x)Il{t,x)dx 

+ — / e^^x^^(t^x)dx. (73) 

fi Jo 

Thus, 


||G(v,^)||, ||0(x)||,||/i(x,^)||<M, 

VvG[0,l],^G[0,x], (69) 

where the matrix/vector norms || • || are compatible with 
the other corresponding matrix/vector norms. Hence, using 
Young’s inequality, the following are derived 


if a 2 e {t,x)rinvOll{t,x) dx 

Jo 

< 2 — f a 2 e~^’‘Yf{t,x)n{t,x)dx 

£ Jo 


(70) 


and. 


fl ^d2X ^ 

/ - ^(t^x) 6 (x)Il(t^x)dx 

Jo M 

/*! 

- / 

Jo M 

+n^(?,x)0^(x)0(x)n(r,x)) dx 

fl g^x _ 

< / - ^^{t,x)dx 

Jo M 

Nf a T 

+ — / (t,x)mt,x)dx 

W Jo 


(71) 


and, 


/ / a 2 e {t,x)TinMx,^) dx 

Jo Jo 

M px 

<— a 2 e~^^ / (n^{t,x)n(t,x) 

£ Jo Jo 

+U'^{t,^)U{t,^))d^dx 

M ^ rj. 

= — / a 2 e (t,x)Il{t,x)dx 

£ Jo 

J a 2 {t,x)Il{t,x)dx 


M fl 
S 2 £ . 


<a 2 e 


!'• 


- 82 X 


M M \ T 

—xH- {t^x)Il{t^x)dx (72) 

£ 02 £ ’ 


V 2 {t)<-e ^a2(7rf(f,l) + #2(^l)) 
+ e^ (Pi7ri(?,l) + P2#2(?,l)f 


7'‘'- 

Jo 


1 X 

'e-^^{& 2 a 2 


^^(t,x) dx 


a 2{2 + x+l/ 52 )M \ _ 

e J &2^^ 

Il{t^x)dx. 




S 2 


M 


With the help of the boundary conditions (55), we obtain 

V2(0 < - [(02 - 2p^2^) 1) 

+ (a2-2p2e^^)7r2(^l) 


/ 


- / €^’‘{52- 


\ X 


^^(t,x) dx 


- [^lf(t,x)e-^^P(x) n(t,x)dx, (74) 

Jo 


where 


P{x) =a2 82 — 


( 2 -\-x-\- 1/ 52 )M \ 


252W l_i \ ^j^_j:^g 52 ii+x)^ 


52 


m 2 m 2 


P 52P 


(75) 


First, from (74), we need to choose the tuning parameter 
^ ^. Then, by choosing 


(22 > max 


{2pre-‘>,2p|e-‘=, 


^ _ (2+.y+l/^2)^ 


(76) 
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to make sure that the matrix P{x),x e [0,1] is positive defi¬ 
nite, we could derive exponential stability of the target error 
system. 


Then, from the continuity and invertibility of the backstep- 
ping transformation (53), we could derive exponential con¬ 
vergence of the designed observer. Thus, the following the¬ 
orem is proved. 

Theorem 2 Under the assumptions that the initial data are 
in (^^([0,1]))^, the observer (48) (with the coefficient func¬ 
tions Pi{x)p = 1,3 determined by (56)-(58) and with the 
boundary condition (49)) exponentially convergent to the 
system (15) in the sense. 


5 Output Feedback Control 


The controller (33) requires a full state measurement and the 
observer is designed to reconstruct the state over the whole 
spatial domain based on an output measurement w(r,0). 
Thus, by combining these two, we could design an observer- 
based output feedback controller. 

Theorem 3 Consider the (u\^ w)^- system (15)-(16) to¬ 
gether with the (wi, w)^- observer (48)-(49). For a given 

initial condition (t/^, u^^ w^, u^^ w^)^ G (^^([0, 1]))^ 

and the control law 


U{t) =-piUi{t,l)-p2U2{t,l) 

Jo L 

-f k2(l,(^)w2(^,^)+^3(l,^)w(l,(^)j (77) 


where ki, k 2 and satisfy (21) with the boundary condition 
(22), the (t/i, U 2 ^ w, wi, U 2 ^ w)^-system is exponentially 
stable in the sense of the .^^-norm. 


PROOF. From the definition of the error variable vector 
(50), the combined closed-loop (wi, U 2 ^ w, wi, ^2, w)^- 
system of (15)-(16), (48)-(49) and (77) is equivalent with 
the (^1, ^2, ^2, w)^-system of (48)-(49), (51)-(52) 

and (77). In comparison to the backstepping transformation 
(17) and (18), the invertible transformation 


^i{t^x) = uft^x) for / = 1, 2 (78) 

X{t,x)=w{t,x)- [ ki{x,^)ui{t,^)d^ 

Jo 

- rk2{x,^)u2{t.^)d^ 

Jo 

- f h{x,^)w{t,^)d^ (79) 

Jo 


and (53) maps the system (48)-(49) into a (i/i, ^ 2 , X, ^1, 
^2, 0)^-system, of which the exponential stability can be 
proved through the following Lyapunov function: 


Jo V ri ji J 

'^-x'^{t,x)dx 


-L 






+ b 


f 


a2e 




ri 


72 


dx 



(80) 


Exponential stability of the (wi, M2, w, mi, U 2 ^ M>)^-system 
is thus proved. 


6 Numerical simulations 

This section is devoted to the numerical simulations of sys¬ 
tem (13) subject to the boundary conditions (16) using re¬ 
spectively the controller U(t) defined in (33) and (79). Our 
goal is to demonstrate the performance of the suggested con¬ 
trollers (33) and (79) to stabilize system (13) around the zero 
equilibrium. For the sake of completeness, we give a short 
description of the used numerical schemes. We employ an 
accurate finite volume scheme to advance in time and space 
the hyperbolic evolutionary system (13). Elsewhere, for the 
implementation of the control law (79), a resolution of the 
kernel PDF’s system (21)-(22) on T is requested. For this 
end, in sight of the triangular shape of T, the finite element 
setups are used. The solution of the kernel problem are com¬ 
puted accurately by using the quadratic finite element pair 

ft. 

The mesh of the triangular domain contains 7655 de¬ 
grees of freedom. For the evolution equation, the com¬ 
putational domain is the segment [0,1] and is divided 
uniformly in 100 cells. Further, the CEL number is fixed 
at 0.9 and the time step 2st in this numerical simulation 
is given by a CEL (Courant-Friedrichs-Lewy) stability 
conditions The initial bottom topography is defined as 

^(0,.x) = 0.4^1+0.25exp with a gaussian 

distribution centered at the middle of the domain. The initial 
water level and its velocity field are computed, respectively 
as Ff(0,v) = 2.5 -5(0,v) and //(0,v)y(0,v) = \{)sin{Kx). 
From the physical variables of F/^(0,v), y(0,v) and 5(0,v), 
the initial data of the characteristic variables v, ui and U 2 
are computed combining relations (6) and (8). It is interest¬ 
ing to mention that these initial conditions imply a strong 
perturbation in the domain. 

6.1 State feedback under subcritical flow regime (ft < 1) 

Let us consider the set point (W, y*, 5*) listed in Table 
1 (see Appendix) and make use of the state feedback con- 
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troller U(t) defined in (33). The chosen set point leads to 
the following characteristic speeds: 

Ai = -1.42, X 2 = 0.76 and A3 = 7.42. 

Besides, the Fronde number is Fr = 0.6 which correspond 
to a subcritical fiow regime. The coefficients at, 6i and the 
matrix a are computed with the help the characteristics 
speeds A/. Hence the kernel PDFs (21)-(22) is solved nu¬ 
merically and the value of the kernel ki , ^2 and at v = 1 
are employed for the implementation of the state feedback 
controller (33). An illustration of the kernel solution ki is 
presented in Figure 4. 



Figure 4. Numerical solution of the kernel component ki on T 


In Figure 5 are depicted the behavior in time of the control 
U{t) and the output measurement y{t). Clearly, despite the 
initial amplitude of U{t), this latter one decreases in time 
and vanishes after t > 4s. Let us remind that the implemen¬ 
tation of U{t) requires a full-state measurement. Moreover, 
likewise U(t) the output measurement y{t) shows the same 
trend with its amplitude decreasing in time and tending to 
zero after t > 3s SiS can be seen in Figure 5(b). 

Plus, in figure (6) we plot the evolution in time of the 
norm of the characteristics. As expected from the theoritical 
part we observe that the norm of the characteristics con¬ 
verge to zero. As a result this shows that the system (15) 
converge to the zero equilibrium. Thereby the physical lin¬ 
earized model (2 ) also converges to (H*, V*, 5*). 

In addition. Figure 7 describes the space and time dynam¬ 
ics of the plant and is consistent with the numerical results 
presented above. As time increases, we notice that the per¬ 
turbation in the overall system decreases and vanishes later. 

6.2 Output feedback under supercritical flow regime (/v > 

1 ) 

Here, all parameters of the physical model are listed in the 
following Table 2 given in the Appendix. In this subsection, 
the dynamic of the closed-loop system (13) together with the 
output feedback control law (79) is simulated The set point 



(a) Output control law 



(b) Measured output 


Figure 5. Evolution in time of the control law U(t) and the mea¬ 
sured output y{t). 



Time 

Figure 6. Evolution in time of the norm of the characteristic so¬ 
lution. 

{W, y*, 5*) leads to the following characteristic velocities 

Ai = 1.87, A2 = -0.74 and A3 = 8.13. 

The Froude number is set to Fr = 1.6. This test case is par¬ 
ticularly challenging since the fiow regime is supercritical. 
As previously we solve the kernel problem As previously 
the kernel problem (21)-(22) is solved numerically and the 
solution is used for the computation of the feedback control 
law (79). Not only that, the system (56)-(57) is also solved 
using the finite element setup and used to compute the kernel 
gain pi{x) defined in (58). This observer gain is represented 
in Figure 9. 
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(a) Evolution of ui{t,x) 
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(b) Evolution of U 2 (t^x) 



(c) Evolution of v(t,x) 


Eigure 7. Behavior in time and space of the characteristic solutions. 



Eigure 8. Solution component k\ on the triangular domain T. 


The value of the kernel ki , and at x = 1 are the gain 
of the designed output feedback controller (79). Elsewhere, 



Eigure 9. Computed observer gains Pi(x). 


In Figure (8) is depicted a snapshot of the numerical solution 
k\ of the kernel PDFs (21)-(22) with the x and coordinates 
defined as the horizontal and the vertical axis, respectively. 


the computation of the control law (79) requires also the 
knowledge of the observer. Then system (54)-(55) is solved 
on time and space. Figure 10 shows the evolution in time of 
the control input U{t) at downstream and the output mea¬ 
surement y{t) at upstream. Clearly, the amplitude of U{t) 
decreases in time and vanishes for t>As. 
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(a) Output control law 



Figure 10. Evolution in time of the control law and the measured 
output 


Moreover, the output measurement y{t) shows the same 
trend with its amplitude decreasing in time and tends to zero 
after t>3. 

The dynamic of the ^^-norm is directly related to the mag¬ 
nitude of the propagation speeds A/ as illustrated in Fig¬ 
ure 11. Furthermore, we give a comparison of our output 
feedback law (Figure 11(a)) to the approach in [11] (Fig¬ 
ure 11(b)) under this supercritical flow regime (fast rapid 
flow) where the conditions of Theorem 2 (cf [11]) are not 
fulfilled. Altogether, our approach exhibits a successful sta¬ 
bilization of system (13) around the zero equilibrium while 
instabilities are noticed when using the strategy presented in 
[11]. Figure 12 describes the space and time dynamics of the 
plant and is consistent with the numerical results presented 
above. As time increases, we notice that the perturbation in 
the overall system decreases and vanishes later. 

As can be seen from these numerical simulations, the sys¬ 
tem (13) subject to the feedback control U{t) is stabilized 
around the zero equilibrium as expected from the theoret¬ 
ical part. We clearly observed the expected qualitative and 
physical behavior of the designed control regardless the na¬ 
ture of the fiow. 

This class of state feedback law and output feedback con¬ 
troller we apply here allows to stabilize the SVE system in 
a optimal time when comparing with the results presented 



Time [s] 


(a) Ouput feedback control through the backstep- 
ping design. 



(b) Lyapunov design when the requirements of 
Theorem 2 in [11] are not fulfilled. 


Figure 11. Evolution in time of the norm of the characteristic 
solution. 

in [12] where a boundary measurement scenario is adopted. 

7 Conclusion and future work 

In this paper, a linearized Saint-Venant Exner model is an¬ 
alyzed for control purposes. The model describes the evo¬ 
lution of the water fiow coupled with the transport of a 
sediment layer in an open channel. A backstepping state 
feedback controller, located at the downstream gate of the 
channel, is first designed for the (exponential) stabilization 
of the water level and the bathymetry at a desired equilib¬ 
rium set, under the subcritical or supercritical fiow regime. 
Then, based on an exponentially convergent Luenberger ob¬ 
server that reconstructs the full state, we design a a back- 
stepping output feedback controller with the measurements 
at upstream. This controller also achieves the exponential 
stability of the linearized SVE model, for both subcritical 
and supercritical fiow regime. Although the backs tepping 
approach offers a more complicated design than the method 
developed in [11], it enables the exponential stabilization of 
the SVE system without any restriction on the system and 
the nature of the fiow. It also reduces the number of actua¬ 
tors of the system: we only need a single boundary control 
in this paper, but on-line measurements of the water levels 
at both ends of the spatial domain are needed in [11]. More- 
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(a) Evolution of ui{t,x) 



(b) Evolution of U 2 {t,x) 



(c) Evolution of v{t,x) 


Eigure 12. Behavior in time and space of the distributed states. 

over, simulation results in comparison to [11] are provided 
to verify that the proposed controller moves beyond those 
limitations of [11]. 

We emphasize that practically, such systems are subjected 
to several types of perturbations and model uncertainties. 
Thus, an effective control action must take into account of 
these factors. We refer the readers to the recent results pro¬ 
posed in [28,27] on the stabilization of hyperbolic PDEs 
with matched disturbances at the boundary input. In these 
papers, the authors employ sliding mode control and active 
disturbance rejection control to deal with them. Our future 
objective is to consider robustness issues for this application. 
Also, the extension of this approach to a network of flow and 
sediment transportation remains an interesting open problem 


with a high potential in real applications. Among the oth¬ 
ers, the problem of proving local stability of the nonlinear 
SVE plant under the linear feedback will be very interested 
to study. 

Appendix 

• Subcritical flow regime state feedback 


T 

Ax 

CFL 

A 

P 

Cf 

Pi 

P2 

8 

0.01 

0.95 

0.008 

0.002 

0.1 

1.5 

1.5 


q\ 

qi 

H* 

u* 

B* 


1 

1.2 

2 

3 

0.4 


Table 1 

Physical parameters and dimensionless numbers 
• Supercritical flow regime output feedback 


T 

Ax 

CFL 


Pg 

Cf 

Pi 

P2 

8 

0.01 

0.9 

0.003 

0.002 

0.1 

1 

1.5 


q\ 

qi 

H* 

u* 

B* 


1 

1.2 

1 

5 

0.4 


Table 2 

Physical parameters and dimensionless numbers 
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